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A large body of experimental evidence suggests that the decay of the false vacuum, accompanied 
by quantum pair creation of soliton domain walls, can occur in a variety of condensed matter 
systems. Examples include nucleation of charge soliton pairs in density waves [eg. J. H. Miller, 
Jr. et al., Phys. Rev. Lett. 84, 1555 (2000)] and flux soliton pairs in long Josephon junctions. 
Recently, Bias and Lemos [J. Math. Phys. 42, 3292 (2001)] have argued that the mass m of the 
soliton should be interpreted as a line density and a surface density, respectively, for (2+l)-D and 
(3+l)-D systems in the expression for the pair production rate. As the transverse dimensions are 
increased and the total mass (energy) becomes large, thermal activation becomes suppressed, so 
quantum processes can dominate even at relatively high temperatures. This paper will discuss both 
experimental evidence and theoretical arguments for the existence of high-temperature collective 
quantum phenomena. 

PACS numbers: 03.75.Lm, 71.45.Lr, 75.30.Fv , 85.25.Cp, 11.27.+d 



I. BACKGROUND AND MOTIVATION 

A wide class of nonperturbative phenomena in field 
theory can be understood in terms of quantum tunneling. 
A well-known example is the quantum decay of the false 
vacuum jwhich has been of broad scientific interest in 
cosmology 0,0 and other fields "3] for over two decades. 
In three dimensions, the boundary between the bubble 
of true vacuum and the surrounding false vacuum is a 
type of topological defect known as a domain wall. A va- 
riety of topological defects in condensed matter systems 
have been proposed to nucleate via quantum or thermal 
fluctuations. These include vortex-antivortex pairs and 
vortex rings in superconductors 0, 0] , superfluids 0, , 
and Bose-Einstein condensates |9|; dislocation pairs in 
Wigner crystals ^^'^ vortex lattices 0; phase -slip 
vortex rings in charge density waves (CDWs) [Ulia; 
charge (or flux) soliton-antisoliton pairs in density waves 
0j O (or Josephson junctions JL^]); and soliton-hke 
domain walls surrounding cigar-shaped bubbles of true 
vacuum in three-dimensional CDWs [TtI ]. 

Topological defects, such as flux vortices, play an es- 
pecially important role in the cuprates and other type-II 
superconductors. Magnetic relaxation rates that depend 
weakly on temperature up to 20 K J^J , or even decrease 
with temperature suggest that vortices may tunnel 
over a wide temperature range. Moreover, the consis- 
tently low IcRn products of cuprate Josephson devices 
suggest that Josephson vortex-antivortex pair creation 
|16| may occur when the current is much smaller than 
the "classical" critical current la ^ A/i?„e. 

A charge density wave (CDW) is a condensate po| 
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in which the electronic charge density in a quasi-one- 
dimensional metal is modulated, p{x,t) — pQ(x,t) + 
piCos[2kpx — (j){x,t)]. Here po{x,t) contains the back- 
ground charge of the condensed electrons, and an excess 
or deficiency of charge proportional to idcjj/dx. A spin 
density wave (SDW) has a modulated spin density, and 
is equivalent to two out-of-phase CDWs for the spin-up 
and spin-down subbands. Although pinned by impuri- 
ties, a density wave (DW) can transport a current when 
an applied field exceeds a threshold value Et- The most 
widely studied models of DW depinning are variations of 
that proposed by Fukuyama, Lee, and Rice pl|. 

A long-standing debate concerns whether a classical 
description is sufficient to describe density wave depin- 
nin g, or a quantum treatment is required. John Bardeen 
|22. 123 ] proposed a model in which condensed, dressed 
electrons Zener tunnel through a "pinning gap," whose 
energy per electron is small compared to the Peierls 
gap. This model was motivated by the Zener-like behav- 
ior, J ^ [E - Et] cxp{-Eo/E) j^l observed in the 
current-field characteristics of NbSea and TaSs. Bardeen 
pointed out that the 3-D coherence of the DW conden- 
sate would suppress thermal excitations without neces- 
sarily suppressing the amplitude for coherent tunneling 
of the internal microscopic degrees of freedom. The ob- 
served Zener field Eq does not increase with the trans- 
verse dimensions, which has important implications for 
any theory based on tunneling probabilities. 

Maki considered a real-space picture, in which 
27r-soliton-antisoliton {S — S') pairs nucleate via quan- 
tum tunneling. Krive and Rozhavsky pointed out 
the existence of a Coulomb blockade threshold field for 
the creation of charged S ~ S' pairs, and later extended 
this picture to nucleation of soliton-like domain walls 
|T3|- A more recent paper [23| proposed an analogy 
to time-correlated single electron tunneling to explain 
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the observed lack of DW polarization below threshold, 
and to interpret key features of DW dynamics, such as 
coherent oscillations and narrow-band noise. Observa- 
tions of Aharonov-Bohm oscillations '2^ in the magneto- 
conductance of CD Ws in NbSc3 crystals with columnar 
defects strongly support the idea that quantum trans- 
port plays a fundamental role. Additional evidence for 
quantum behavior in density waves is found in rf ex- 
periments [2^ that show very good agreement with 
photon-assisted tunneling theory at temperatures up to 
215 K. 

It is important to also consider whether or not den- 
sity wave dynamics can be interpreted classically. Per- 
haps the most obvious requirement for classical sliding 
is that the washboard potential should be tilted enough 
for the deformable object to slide over the pinning bar- 
rier. This tilting of the washboard potential leads to a 
predicted phase displacement of 7r/2 in the s-G model 
at the classical threshold, and to even greater predicted 
phase displacements when disorder is included jsj . How- 
ever, NMR and high-resolution x-ray scattering experi- 
ments show that the phase displacements of both charge 
[3^ I33L Is^ l and spin [35| density waves are more than 
an order of magnitude smaller than classically predicted 
values near threshold. 

Classical models also predict |^ |3^ that the low- 
frequency dielectric response should increase, or even di- 
verge, as the classical threshold is approached from be- 
low. Such predictions are refuted by the observed bias- 
independent rf and microwave responses below threshold 
in both CDWs US 113 and Wigner crystals These 
experiments strongly suggest that the washboard pinning 
potential is barely tilted as the measured threshold field 
is attained. Thus, the observed threshold appears to be 
much smaller than the classical depinning field, and may 
represent a Coulomb blockade threshold for the creation 
of charged topological defects. Moreover, attempts [s^ 
using a scanning tunneling microscope (STM) to directly 
observe either displacement of the CDW below thresh- 
old or sliding above threshold have been unsuccessful. 
The apparent lack of sliding seen in STM experiments 
and the jerky dynamics revealed by NMR experiments 
[3^ I33I I suggest that the DW spends most of its time in 
the pinned state even above threshold. 

The classical deformable DW model also predicts that 
the high-field dc conductivity should have the form 
o'cdw ^ a — hE^^/"^ , where a and b are constants This 
prediction is contradicted by experiments that show 
substantial departure from the classically predicted be- 
havior, but essentially perfect agreement with the Zener 
form ^ exp(— i?o/i?), in the high-field limit. Copper- 
smith and Littlcwood '4^ (CL) make the following pre- 
dictions for mode locking of a moving DW to an ac source 
within the classical model. (1) Damped relaxation is cru- 
cial in that complete mode locking can occur only when 
the ac frequency and amplitude are small enough for sig- 
nificant relaxation to take place when the total field is 
less than Et- (2) For large ac amplitudes and frequen- 



cies such that the time interval spent below threshold 
is small, the differential resistance dV/dl should exhibit 
peaks but not complete mode locking. In addition, no 
"wings" (sharp negative dips in dV/dl adjacent to the 
interference peak) should be observed. These predictions 
are refuted incontrovertibly by ac-dc interference and 
mode locking experiments on high quality NbSca crystals 
41] that demonstrate both negative wings and complete 
mode locking at high dc bias fields and high frequencies. 

The failure of CLs calculations to account for the 
mode-locking experiments reflects a fundamental defi- 
ciency of the classical deformable model. When the DW 
slides rapidly, relaxation into the pinning potential wells 
cannot occur quickly enough to prevent the deformations 
and net pinning energy from vanishing. The character- 
istic frequency for this relaxation, the DW dielectric re- 
laxation frequency, is measured to be about 5 MHz in 
NbSca 01 . The observations of complete mode locking 
at ac and drift frequencies of up to 1 GHz indicate that 
the strength of DW pinning is undiminished in a regime 
where significant relaxation cannot occur within the con- 
ventional classical model. Thus, quantum mechanisms 
must be explored to understand, fundamentally, DW de- 
pinning and dynamics. 



II. SOLITON TUNNELING MODEL 

A realistic density wave pinning potential would have 
the form J7o(r){l — cos[(^ — 0o(r)]}, which includes spatial 
variations of both the pinning energy J7o, and the opti- 
mum phase (j)Q. However, observations of narrow band 
noise, coherent oscillations, and complete mode-locking 
with an ac source 41'] in high quality NbSea crystals 
suggest that both [/q a-nd 0o are slowly varying in weakly 
pinned DWs. The sine-Gordon (s-G) model thus provides 
an idealized conceptual frameworkjfor interpreting many 
aspects of density wave transport |4j| , from both classical 
and quantum points of view. Similarly, a weakly-coupled 
Josephson junction (JJ) can be described by a phase, 
(j){x, t), representing the phase difference across the junc- 
tion, in a sine-Gordon potential. 

Density waves and JJs are dual in that the roles of 
charge and flux are interchanged, as well as those of 
current and voltage. The current in a density wave is 
/ = {Qo/2TT)d(l)/dt, where Qq ~ 2eNch and Nch is the 
number of parallel chains, whereas the voltage across a 
JJ is y = {<^o/'2Tr)d(j)/dt, where $0 = h/2e. Charge 
(flux) solitons in a density wave (JJ) carry a charge (flux) 
of ±Qo{±^q) (see Table HJ. The width of a Joseph- 
son vortex is roughly the Josephson penetration length, 

— 1/2 

Xj OH Jc whereas, in a DW, the soliton width is 
'^o — cq/^Oi where cq is the phason velocity and ujq is 
the pinning frequency. Thus, Aq will increase with de- 
creasing impurity concentration (as loq decreases), and 
may approach the distance between contacts in extremely 
pure samples. This is equivalent to approaching the short 
junction limit, L < Aj, in a JJ, where the V — I curves 
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TABLE L Charge-flux duality between density waves and Josephson junctions. 



Density Wave Josephson junction 

Sohton or antisoliton Kink w/ charge ±Qo Josephson vortex w/ flux ±<^o 

Type of threshold Threshold fleld Et Threshold current It 

Transport characteristic I vs. V , I — 2TiQod<j)/ dt V vs. I, V — 2Ti(f)iid(j>/dt 



become significaiitly less rounded. 

The quantum decay of the metastable false vacuum 
of a scalar field (/), accompanied by the creation of soli- 
tons and antisolitons in (1+1) dimensions [islli^li^ has 
been studied extensively in quantum field theory. Dias 
and Lemos (DL) have calculated the effective one 
loop action for pair creation of solitons in the (l-l-l)-D 
s-G model. In addition, they have extended their argu- 
ments to pair creation of soliton domain walls in higher 
dimensions. In the (l-l-l)-D case, the soliton-antisoliton 
pair production rate per unit time and length, F/i, is 
identified with the decay rate of the false vacuum. This 
is calculated to be, in the absence of Coulomb interac- 
tions between the solitons (taking h ~ c ^ 1) [isf: 



- = — exp — — (Ij 
L In L s J 

where m is the soliton mass (energy). Here e represents 
the gain in energy per unit length between the two cre- 
ated solitons, which is proportional to the applied electric 
field E for a density wave (and to the applied current for 
a JJ). The energy for pair creation is contained in the line 
between the particles and is given by 2xs = 2m, where 
2x is the distance between the solitons. 

DL point out ji^l that a one-particle system in (1+1 )-D 
can be transformed into a line in (2+l)-D and a thin wall 
in (3+l)-D, except that the mass m should be interpreted 
as a line density and a surface density, respectively. Their 
calculations thus apply directly to the domain wall pair 
creation. As the transverse dimensions are increased and 
the total mass (energy) becomes large, thermal activation 
becomes suppressed, so quantum processes can dominate 
even at relatively high temperatures (up to 400 K for the 
CDW in NbSea). 

As an example, consider the (2+l)-D case, where one 
can make a cut normal to the infinite strings so that "e" 
will still be the energy density of the line between the 
two infinitesimal pieces of the strings. Each one of these 
infinitesimal pieces has a mass m, so m represents the 
line density of the string. The energy for pair creation 
of strings comes from the infinite plane between the two 
strings, and is given by 2xde — d2m, where d is the length 
of the string, which goes to infinity jl^ . Therefore both 
the mass of the string and the energy to create the pair 
of strings become infinite. 

Importantly, "m" can no longer be interpreted as the 



total mass of the string or of the domain wall. Indeed 
such an interpretation would correspond to creation of 
point particles in (2+l)-D or in (3+l)-D, where the total 
mass m would be infinite. The infinite energy for pair 
creation would then have to come from the line between 
the two particles rather than from the infinite plane be- 
tween the two strings. DLs calculations do not apply to 
this process. In the generalization to higher dimensions, 
one is using the symmetry of the system along the ex- 
tra direction. In Minkowski space-time, the procedure 
is direct However, in systems that include gravita- 
tion, the procedure is usually not so straightforward, and 
one must employ dimensional reduction or dimensional 
expansion jSOj. 

Suppose we now put the constants h and c back into 
Eq. ^and take A to represent the Compton wavelength, 
A = h/mc. (In a density wave, the effective "speed of 
light" is the phason velocity co, which is about 3 x 10^ 
cm/s.) Explicitly putting the physical constants into the 
exponent, we then have: rmr'^/e — > rmrc^ /[qEX], where 
q is the soliton charge for a single chain (i.e. for the 
(1+1)-D system) and E is the electric field. (In a JJ, 
the charge and electric field would be replaced by the 
fiux $0 and the current density.) For the (3+l)-D sys- 
tem, we can now define m' = m/a^ and q' — q/a^ to 
be the mass and charge per unit area of the soliton do- 
main wall, respectively. In a density wave, the length 
represents the transverse distance between DW chains. 
Note that, for the exponent Tim' /[q' EX\ to remain di- 
mensionless, A needs to retain its original units of length. 
One way of doing this is to define a new "quantum of ac- 
tion per unit area" in terms of the characteristic length 
scale a, i.e. h! = hja} , leading to a Compton wave- 
length, A = h! /m'c, which has the same value as before. 
Lattice gauge theorists have recently examined models of 
domain wall fermions (bosons), whose anticommutators 
(commutators) differ from those of point particles, thus 
lending credence to the arguments presented here. 

Another way of looking at the problem is to recall that 
h has units of angular momentum: L — rmv, where r 
= distance, m — mass, and v — velocity. Now, since 
the mass transforms as m' — m/A, where A is an area, 
from L = rmv we must have U — L/A, so that h' — 
h/A The commutator between bosonic domain wall 
creation and annihilation operators will then be scaled 
accordingly. The Compton wavelength of the (3+l)-D 
problem refers to the direction normal to the domain 



wall, and is thus A = h! /m'c = h/mc. 

It should be noted that rf experiments 
are consistent with a ratio h/q = h! /q' that is scale in- 
variant, where q = 2e for a fully condensed CDW. This 
is further supported by magneto-transport experiments 
[28| , which yield Aharonov-Bohm oscillations with a peri- 
odicity of h/2e, and not h/2Ne as predicted theoretically 
|5l|. where N is the number of coupled chains. However, 
the observed mode locking with an ac source 141]] shows 
that the phase is coherent throughout macroscopic re- 
gions within the crystal, suggesting quantum nucleation 
of entire domain walls rather than just single dislocations. 

The above arguments lead, for a CDW, to a Zener-like 
current-field characteristic given by: 




J Ee^pi-Eo/E] 



(2) 



where the characteristic Zener field Eq = [Trm''^ Cq] / [h! q'] 
is independent of the traverse dimensions. (For a JJ, 
the V vs. I curve would be given by y ~ / exp[— /q/^]-) 
DLs calculations neglect Coulomb interactions between 
charge solitons (or magnetic interactions between flux 
solitons in a J J) , and thus fail to predict a sharp thresh- 
old field (threshold current for a JJ) for pair creation. 
The following section examines the origin of the thresh- 
old field as a macroscopic Coulomb blockade effect (or its 
dual for the case of a JJ). 



III. MACROSCOPIC COULOMB BLOCKADE 
AND TIME-CORRELATED SOLITON 
TUNNELING 

The quantum interpretation of the threshold field, as a 
pair-creation threshold due to Coulomb blockade, is mo- 
tivated by Colemans paper [5^ on soliton pair-creation 
in the massive Schwinger model. A pair of S and S' do- 
main walls with charges iQo produce an internal field 
of magnitude E* — Qq/eA, as shown in Fig. ^ where 
A is the cross-sectional area and e is the dielectric con- 
stant. When a field E is applied, the difference in electro- 
static energies of a state with a pair of separation / and 
of the "vacuum" is AU = {l/2)e Al[{E + E*f - E"^] = 
Qol[{l/2)E*±E] , which is positive when |£;| < {l/2)E*. 
Conservation of energy thus forbids pair production for 
fields less than a quantum threshold, Et = (l/2)i?* — 
Qo/2e^, which appears to be about two orders of magni- 
tude smaller than the classical depinning field in NbSca 
[23. The observed universality relation, eEt ~ eNch/A 
|20l| . thus arises quite elegantly in this model. Screen- 
ing by the uncondensed normal carriers greatly enhances 
the dielectric response of NbSca , which has an incomplete 
Peierls gap and also has a much lower threshold field than 
most other CDW materials. (In a JJ, charge-fiux dual- 
ity leads to a quantum threshold current proportional 
to <^q/2Lm, where Lm is the mutual inductive coupling 
between the Josephson vortex and the antivortex.) 



FIG. 1: (a) Density wave phase vs. position, illustrating 
the production of a soliton-antisoliton domain wall pair, (b) 
Model of a density wave capacitor, showing the nucleated do- 
main walls moving towards the contacts. The electric field 
between the domain walls is reduced by the internal field E* . 
The distances Z, Ao, and the crystal thickness are greatly ex- 
aggerated for clarity. 



A density wave between two contacts behaves as a ca- 
pacitor with an enormous dielectric constant, as shown 
in Fig. n] The initial charging energy is (5^/2C, where 
Q is the displacement charge and C = eA/L. We define 
e = 27rQ/Qo = 2t: E/E* = ttE/Et and note that a dis- 
placement (j) near the middle creates a non-topological 
kink-antikink pair, with charges i((j)/2n)QQ, if (j) — 
at the contacts. The washboard pinning and quadratic 
charging energies can then be written as |27] : 



U[4>] 



dx {up[l — cos (/'(x)] + Uc[9 — (j){x)Y} (3) 



where the first and last terms represent the pinning and 
electrostatic charging energies, respectively, and where 
Up ^ Uc for NbSca- If the system starts out in its ground 
state, conservation of energy will prevent tunneling when 
the applied field is below threshold, 9 < tt {E < Et), as 
illustrated in Fig. |21 However, when 6 exceeds tt, what 
was formerly the true vacuum becomes the unstable false 
vacuum. One or more bubbles of true vacuum, with soli- 
ton domain walls at their surfaces, then nucleate and 
expand rapidly (Fig. QJ. After n solitons of charge Qq 
(and antisolitons of charge —Qq) have reached the con- 
tacts, the charging energy becomes: 



(Q-nQo)' 
2C 



27rn)2 



(4) 



This series of piecewise parabolas is similar to the 
charging energy of a single-electron tunnel junction, ex- 
cept that Qo now represents a macroscopic charge. 
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FIG. 2: Plot of potential energy vs. <j) for two different values 
of 6, in which the many spatial degrees of freedom are illus- 
trated schematically. Tunneling is prevented by conservation 
of energy when 9 < n. When 9 > tv, parts of the system 
tunnel into the adjacent well via decay of the false vacuum. 



The single-electron transistor (SET) ^] consists of 
a gate capacitor Cg coupled to an island electrode be- 
tween two small capacitance tunnel junctions in series. 
The gate voltage modulates the I — V curves between 
the source and drain electrodes, with a period e in 
displacement charge, Qg — CgVg. The displacement 



charges Qi,2 across the two tunnel junctions are related 
as Q2 — Qi + Qg + Qo, where go is a phenomenologi- 
cal offset charge induced during cooling The SET 

is related by charge-flux duality to the dc SQUID. The 
critical voltage across an SET is a periodic function of 
Qg, whereas the critical current across a SQUID is peri- 
odically modulated (with period <i>o) by the flux $. 

The model discussed above suggests that it may be 
possible to demonstrate a macroscopic version of the 
SET, by attaching a gate capacitor to an island electrode 
near the center of a quasi- 1-D crystal with a density wave. 
The displacement charge induced by the gate electrode 
would then periodically modulate the total critical volt- 
age between the source and drain electrodes. Ideally, in 
the absence of any shunt conductance, the periodicity of 
the gate displacement charge might be expected to be 
~ Qq. However, screening by the normal, uncondensed 
electrons will tend to reduce the effectiveness of the gate 
which, unlike the source and drain contacts, cannot be 
driven by a current source. The displacement charges 
across the two segments of the crystal will be related as 
Q2 = Qi + PiQg + 9o), where <C 1 reflects screening 
by the normal carriers. The total charging energy of the 
two segments, in this idealized model, will then be: 



(Qi-niQo)^ , ((32-'^2Qo) 



2Ci 



2C2 



r 



2Ci 



2C2 



(5) 



where Ci and C2 are the capacitances of the two seg- 
ments separated by the island electrode. The analogy 
to the SET suggests that a gate voltage might modu- 
late the I — V curves between source and drain contacts, 
with a periodicity AVg ^ Qo/PCg- The gate capaci- 
tance Cg, and hence the attainable displacement charge 
Qg, may have been too small to observe non-monotonic 
behavior in previous experiments |54j . in which a gate 
electrode was fabricated directly on the crystal to form a 
MOSFET-like structure. The soliton tunneling transis- 
tor (STT), discussed in the next section, employs a much 
larger, 1 /iF, gate capacitor coupled to an NbSea crystal, 
and exhibits non-monotonic behavior. 

An exact calculation of the charging energy would yield 
a plot in which the energy is reduced slightly at the cross- 
ing points, 9 = niT, as compared to piecewise parabolas. 
However, Eq. ^ provides a reasonable approximation to 



expanded as a Fourier series: 



E{9) - (61 - 27m)2 



(6) 



As in the SET, the voltage across each segment of 
the STT is related to the charging energy as Vi_2 = 
dE/dQi^2 K dE/d9i^2- If we use the approximation given 
by Eq. |S1 this yields a sawtooth function, which can be 



Vf) saw( 



'1.2) = 



^ sin(n0i,2) 



(7) 



This model can also been used to model the dynamics 
of density waves and to interpret the narrow-band noise 
spectra, the many harmonics of which arc consistent with 
a sawtooth function. 

Equation |7| is related by charge-flux duality to the 
current-phase relation of a Josephson junction. The pe- 
riodic behavior of the source-to-drain voltage vs. gate 
voltage of an STT can readily be understood by exploit- 
ing the duality with a dc SQUID. When each J J in a 
dc SQUID has an ideal, sinusoidal current-phase relation 
and the critical currents Iq are identical, the total current 
is / = /o[sin(/3i -|- sin(/32], where (^2 = + 2<i>/$o- This 
yields a total critical current, Ic = 2/o| cos(27r$/$o)|j 
which is a periodic function of the flux $. Similarly, 
the total source-to-drain voltage of an ideal STT, V = 
V\ + V2, where Vi^2 are given by Eq. [T] yields a crit- 
ical voltage that is a periodic function of gate voltage 
Vg — Qg/Cg, 3.8 shown lu Fig. 13 Here we note that 
92 = 9i + 9g + 9o, where Og = 27r/3Qg/Qo = 27r/3go/Qo, 
Qg is the displacement charge across the gate capacitor, 
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FIG. 3: Predicted critical voltage vs. normalized gate voltage 
6g for several values of using the idealized model discussed 
in Section III, showing the periodic behavior. 



go is the offset charge, and [3 is the screening parameter 
due to the normal electrons. 



IV. SOLITON TUNNELING TRANSISTOR 
EXPERIMENT 

The experiment reported here employs a geometry 
analogous to that of the SET, in which the "source-to- 
drain" I — V characteristic is modulated by a voltage ap- 
plied to a gate capacitor. Single crystals of NbSea were 
employed in the experiment. This material forms two in- 
dependent CDWs, at Peierls transition temperatures of 
145 K and 59 K |2J|, respectively. The Peierls gap opens 
up over most of the Fermi surface (FS) below the lower 
transition, but leaves a small portion of the FS intact, 
so that a significant concentration (~ 6 x 10~"'^*cm~^) of 
normal, uncondensed carriers remain down to low tem- 
peratures. 

The geometry used in our experiment is illustrated in 
the inset to Fig. ^ where the width of the crystal is 
exaggerated for clarity. The NbSea crystal was placed 
onto an alumina substrate with a series of evaporated, 
25 ^m wide gold contacts. The substrate was thermally 
anchored to a cold-finger in the vacuum shroud of an 
open cycle helium flow cryostat and the temperature was 
controlled using a Lake Shore temperature controller at- 
tached to a heater coil wrapped around the cold-finger. 
A Keithley programmable dc current source injected the 
current between two contacts, which were bonded to the 
crystal near the ends using silver paint. The "source-to- 
drain" voltage was measured between two additional gold 
contacts, as illustrated, and a 1 /xF gate capacitor was at- 
tached to the center gold contact using silver paint. The 
spacing between contacts along the crystal was 500 /xm 
center-to-center, and the gate capacitor was kept inside 
the cryostat. 
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FIG. 4: CDW current vs. source-to-drain voltage in a soliton 
tunneling transistor (NbSea) for several values of gate voltage 
at 35K. (The shunt current of the normal electrons has been 
subtracted for clarity.) 



We found that substantially smaller gate capacitors 
(as well as gate capacitors with longer leads) were un- 
able to induce a periodic modulation of the I — V char- 
acteristic. Moreover, dc I — V (rather than differential 
dV/dl) measurements were necessary to avoid inducing 
a displacement current through the gate capacitor. A 
programmable voltage source was coupled to the gate 
capacitor via a lO-kfi resistor, which limited the current 
flowing through the crystal during changes in gate volt- 
age when the gate capacitor either partially charged or 
discharged. The cryostat was kept inside an electromag- 
netically shielded enclosure, and Vsd was measured with 
a nanovoltmeter. 

The measurements were primarily carried out at 35 K. 
Previous "field effect transistor" experiments showed 
the greatest modulation at around 30 K, where the 
threshold field is near its minimum. We also observed 
the largest modulation, using our geometry, at compa- 
rable temperatures. We attained the best temperature 
stability (better than ±0.01 K) when the temperature 
was set to 35 K, and thus chose this temperature for 
most of the measurements reported here. 

Figure0]shows several plots of CDW current as a func- 
tion of source-to-drain voltage, Icdw vs. Vsd, in a NbSes 
crystal at 35 K, for different values of gate voltage Vg. 
The gate voltage is seen to modulate the threshold volt- 
age in the I — V curves of Fig. 21 Figure \E\ displays 
plots of source-to-drain voltage Vsd vs. Vg for three val- 
ues of total bias current above threshold. The plots ex- 
hibit roughly periodic behavior, similar to that observed 
in SETs. However, in our system, the measured period- 
icity AV^ ^ lOV is consistent with a macroscopic dis- 
placement charge AQ = CgAVg ~ 6 x 10^'^e, comparable 
to the charge of the conducting electrons between the 
contacts. 

The behavior shown in Fig. [S] is quite extraordinary. 
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FIG. 5: Source-to-drain voltage Vsd vs. Vg for fixed values of 
total bias current at 35 K. 

and appears to be consistent with the soliton domain 
wall tunneling hypothesis. The number of parallel CDW 
chains, Nch, is about 10^. Thus, one might estimate 
the screening parameter as follows: (3 ~ Q0//S.Q ~ 
2Nche-/^Q 3 X 10~^. However, an alternative in- 
terpretation for the observed periodicity might be that 
all of the normal electrons between the contacts partici- 
pate in screening out the displacement charge Qq of the 
CDW. Thus, the observation that AQ/e ~ 6 x 10^^ is 
roughly the number of conducting electrons between con- 
tacts may not be a coincidence. Further work is needed 
to better understand the effects of screening by the nor- 
mal, uncondensed electrons. 

V. IMPLICATIONS 

The implications of the model and experimental re- 
sults reported here are potentially far-reaching, and could 



impact fields as diverse as cosmology, condensed matter 
physics, quantum computation, and biophysics. In the 
field of cosmology, quantum nucleation of dilaton black 
hole pairs [s^, cosmic strings [H^, and even the entire 
universe 0,13 have been proposed. The European Sci- 
ence Foundation has recently funded a program, known 
as Cosmology in the Laboratory (COSLAB), intended to 
explore laboratory analogs of cosmological objects. Top- 
ics being studied include analogies between vortices in su- 
perconductors and cosmic strings, and artificial (acoustic 
or optical) black holes. 

In the area of biophysics, Frohlich suggested js^, 113 
that long-range quantum coherence may play a funda- 
mental role in biological systems. This may provide a 
mechanism to integrate complex interactions that take 
place within a live cell. A recent model (53|] proposes that 
microtubules, which exist in the cells of all higher organ- 
isms but are especially concentrated in neurons, engage in 
a form of quantum computation. Roger Penrose ,60] has 
even made compelling arguments that large-scale quan- 
tum coherence, perhaps mediated by the microtubules, 
plays a fundamental role in consciousness. It is clear 
that, for Penroses arguments to be valid, some mecha- 
nism is needed to suppress decoherence. Perhaps soliton 
domain walls, or related topological defects, may prove 
to be topologically protected from decoherence at biolog- 
ical temperatures. Further research is clearly warranted 
to explore this extraordinary possibility. 



Acknowledgments 



The authors gratefully acknowledge the assistance of 
James Claycomb, L.-M. (Patrick) Xie, and John Mc- 
Carten as well as valuable insights provided by Oscar 
Dias and Jose Lemos. This work was supported, in 
part, by the State of Texas through the Texas Center 
for Superconductivity and Advanced Materials and the 
Texas Higher Education Coordinating Board Advanced 
Research Program (ARP), and by the R. A. Welch Foun- 
dation (E-1221). 



[1] S. Coleman, Phys. Rev. D 15, 2929 (1977). 

[2] E. P. Tryon, Nature, 246, 396 (1973). 

[3] A. D. Linde, Sov. Phys. JETP 60, 211-213 (1984); Zh. 

Eksp. Teor. Fiz. 87, 369-374 (1984). 
[4] A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 

211 (1981). 

[5] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973). 



[6] Y. Duan, H. Zhang, and S. Li, Phys. Rev. B 58, 125 
(1998). 

[7] G. W. Rayfield, and F. Reif, Phys. Rev. 136, A 1194 
(1964). 

[8] G. E. Volovik, JETP Lett. 15, 81 (1972). 

[9] B. P. Anderson, P. C. Haljan, C.A. Regal, D. L., Feder, 

C. W., Collins, and E. A. Cornell, Phys. Rev. Lett. 86, 

2926 (2001). 



8 



S. T. Chui, Phys. Lett. A 180, 149 (1993). 

G. Blatter, M. V. Feigelman, V. B. Geshkenbein, A. I. 
Larkin, and V. M Vinokur, Rev. Mod. Phys. 66, 1125 [35 
(1994). 

S. V. Zaitsev-Zotov, Phys. Rev. Lett. 71, 605 (1993); S. [36 

V. Zaitsev-Zotov, G. Remenyi, and P. Monceau, Phys. 

Rev. B 56, 6388 (1997). [37 

H. Matsukawa, H. Miyakc, M. Yumoto, and H. [38 
Fukuyama, J. Phys. IV France 9, PrlO-161 (1999). 

K. Maki, Phys. Rev. Lett. 39, 46 (1977); Phys. Rev. B [39 
18, 1641 (1978). 

A. Maiti and J. H. Miller, Jr., Phys. Rev. B 43 12205 
(1991). [40 
A. Widom and Y. Srivastava, Phys. Lett. A 114, 337 
(1986). [41 

I. V. Krive and A. S. Rozhavskii, Soviet Physics JETP 

69, 552 (1989). [42 
H. H. Wen, R. T. Lu, and S. L. Yan (unpublished). 
Y. Y. Xue, Z. J. Huang, P. H. Hor, and C. W. Chu, Phys. [43 
Rev. B 43, 13598 (1991). 

G. Griiner, Density Waves m SoUds, Addison- Wesley, [44 
Reading, Mass. (1994). 

H. Fukuyama, J. Phys. Soc. Jpn. 41, 513 (1976); H. [45 
Fukuyama and P. A. Lee, Phys. Rev. B 17, 535 (1978); [46^ 
P. A. Lee and T. M. Rice, ibid. B 19, 3970 (1979). [47 
John Bardecn, Phys. Rev. Lett. 42, 1498 (1979); ibid. [48 
45, 1978 (1979). 

John Bardeen, Phys. Rev. Lett. 55, 1010 (1985). [49 

P. Monceau, N. P. Oug, A. M. Portis, A. Meersehaut, [50 

and J. Rouxel, Phys. Rev. Lett. 37, 602 (1976); N. P. 

Ong and P. Monceau, Phys. Rev. B 16, 3443 (1977). [51 

R. E. Thorne, J. H. Miller, Jr., W. G. Lyons, J. W. Lyd- 

ing, and J. R. Tucker, Phys. Rev. Lett. 55, 1006 (1985). [52 

I. V. Krive and A. S. Rozhavsky, Solid State Commun. [53 
55, 691 (1985). [54 
J. H. Miller, Jr., C. Ordonez, and E. Prodan, Phys. Rev. 

Lett. 84, 1555 (2000). [55 
Yu L Latyshev, O. Laborde, P. Monceau, and S. 
Klaumnzer, Phys. Rev. Lett. 78, 919 (1997). [56 
J. H. Miller, Jr., J. Richard, J. R. Tucker, and John 
Bardeen, Phys. Rev. Lett. 51, 1592 (1983). [57 
J. H. MiUer, Jr., R. E. Thorne, W. G. Lyons, J. R. Tucker, [58^ 
and John Bardeen, Phys. Rev. B 31, 5229 (1985). [59 
P. B. Littlewood, Phys. Rev. B 33, 6694 (1986). 
J. H. Ross, Jr., Z. Wang, and C. P. Shchter, Phys. Rev. [60] 
Lett. 56, 663-666 (1986); 

J. H. Ross, Jr., Z. Wang, and C. P. Slichter, Phys. Rev. 
B 41, 2722 (1990). 

H. Requardt, F. Ya. Nad, P. Monceau, R. Currat, J. E. 



Lorenzo, S. Brazovskii, N. Kirova, G. Grbel, and Ch. 

Vettier, Phys. Rev. Lett. 80, 5631 (1998). 

W. H. Wong, M. E. Hanson, W. G. Clark, B. Alavi, and 

G. Grner, Phys. Rev. Lett. 72, 2640 (1994). 

O. Narayan and A. A. Middleton, Phys. Rev. B 49, 244 

(1994). 

A. Zettl, and G. Griiner, Phys. Rev. B 29, 755 (1984). 
C.-C. Li, L. W. Engel, J. Yoon, D. C. Tsui, and M. 
Shayegan, J. Phys. IV France 9, PrlO-215 (1999). 

G. Gammie, J. S. Hubacek, S. L. Skala, R. T. Brocken- 
brough, J. R. Tucker, and J. W. Lyding, Phys. Rev. B 
40, 11965 (1989). 

L. Sneddon, M. C. Cross, and D. S. Fisher, Phys. Rev. 
Lett. 49, 292 (1982). 

R. E. Thorne, J. R. Tucker, and John Bardeen, Phys. 
Rev. Lett. 58, 828 (1987). 

S. N. Coppersmith and P. B. Littlewood, Phys. Rev. Lett. 
57, 1927 (1986). 

J. H. Miller, Jr., J. Richard, R. E. Thorne, W. G. Lyons, 

and J. R. Tucker, Phys. Rev. B 29, 2328 (1984). 

M. J. Rice, A. R. Bishop, J. A. Krumhansl, and S. E. 

TruUinger, Phys. Rev. Lett. 58, 828 (1987). 

M. Stone, Phys. Lett. B 67, 186 (1977). 

M. B. Voloshin, Sov. J. Nucl. Phys. 42, 644 (1985). 

V. G. Kiselev, Phys. Rev. D 45, 2929 (1992). 

O. J. C. Dias and J. P. S. Lemos, J. Math. Phys. 42, 

3292 (2001). 

O. J. C. Dias and J. P. S. Lemos, private communication. 
J. P. S. Lemos and V. T. Zanchin, Phys. Rev. D 53, 4684 
(1996). 

E. N. Bogachek, I. V. Krive, I. O. Kulik, and A. S. 

Rozhavsky, Phys. Rev. B 42, 7614 (1990). 

S. Coleman, Annals of Phys. 101, 239 (1976). 

M. H. Devoret, and R. J., Nature 406, 1039 (2000). 

T. L. Adelman, S. V. Zaitsev-Zotov, and R. E. Thorne, 

Phys. Rev. Lett. 74, 5264 (1995). 

F. Dowker, J. P. Gauntlett, D. A. Kastor, and J. 
Traschen, Phys. Rev. D 49, 2909 (1994). 

R. Basu, A. H. Guth, and A. Vilenkin, Phys. Rev. D 44, 
340 (1991). 

H. Frohlich, Nature 228, 1093 (1970). 

H. Frohlich, Proc. Natl. Acad. Sci. USA 72, 4211 (1975). 
N. E. Mavromatos and D. V. Nanopoulos, Int. J. Mod. 
Phys. B 12, 517 (1998). 

Roger Penrose, Shadows of the mind: a search for the 
missing science of consciousness, (Oxford University 
Press, Oxford, 1994). 



